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Abstract. In this paper we are concerned with three lattice problems: the lattice packing 
problem, the lattice covering problem and the lattice packing-covering problem. One way 
to find optimal lattices for these problems is to enumerate all finitely many, locally optimal 
lattices. For the lattice packing problem there are two classical algorithms going back 
to Minkowski and Voronoi. For the covering and for the packing-covering problem we 
propose new algorithms. 

Here we give a brief survey about these approaches. We report on some recent com- 
puter based computations where we were able to reproduce and partially extend the known 
classification of locally optimal lattices. Furthermore we found new record breaking cov- 
ering and packing-covering lattices. We describe several methods with examples to show 
that a lattice is a locally optimal solution to one of the three problems. 



1. Introduction 

Classical problems in the geometry of numbers are the determination of most econom- 
ical lattice sphere packings and coverings of the Euclidean d-space R'^. A lattice i is a 
full rank, discrete subgroup of R"*. Thus there exist matrices A G GLd{M.) with L = ATU^ 
which we call bases of L. 

If B'^ denotes the Euclidean unit ball, then the Minkowski sum L + aB'^ ~ \v ^ ax : 
V G L,x £ B'^}, a E K>o, is a lattice packing if the translates of aB"^ have mutually 
disjoint interiors and a lattice covering if M'' = L + aB'^. The packing radius X{L) of a 
lattice L is given by 

A(i) = max{A : L + XB'^ is a lattice packing}, 

and the covering radius ^{L) by 

~ min{fj, ; L + jiB'^ is a lattice covering}. 

For a lattice L we define its determinant det(i) — \ det(A)|, which is independent of 
the chosen basis. We consider the following three "quality measures" of L: 

(1) the packing density S(L) — — - • k^, 

det[L) 

(2) the covering density Q(L) — — - • k^, 

aet(L) 

(3) the packing-covering constant "f[L) = — — -. 

X[L) 

Here — n'^^'^ /T{d/2 + 1) denotes the volume of the unit ball B'^. So the packing 
density S{L) for instance gives the ratio of space covered by spheres in the lattice packing 
L + X{L)B'^. Note that all three quantities are invariant with respect to a scaling aL of L 
with a 7^ 0. For each of the quantities we consider the problem of finding extremal lattices 
attaining a maximum or minimum respectively. 
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Problem 1.1 (Lattice Packing Problem). For d > 2, determine Sij, — max 6{L) and lattices 
L attaining it. 

Problem 1.2 (Lattice Covering Problem). For d > 2, determine = minO(L) and 
lattices L attaining it. 

Problem 1.3 (Lattice Packing-Covering Problem). For d > 2, determine — min7(L) 
and lattices L attaining it. 

Note that all three optima are attained. All three problems have in common that there 
exist only finitely many local optima for every d (see Section |3] for definitions). In this 
article we want to review some of the major tools available to find such local extrema and 
to verify their local optimality. In Section |2] we briefly summarize known results of the 
three problems and in Section |3] we give a short introduction to the connection of lattices 
and positive definite quadratic forms, which gives the framework in that the problems are 
usually dealt with. In Section |4] we describe two classical approaches by Minkowski and 
Voronoi to enumerate all local optima of the lattice packing problem. In Section |5] we 
are concerned with local optima of the lattice covering and the lattice packing-covering 
problem, which can be treated in parallel. Most of these techniques are described in greater 
detail in l,SV04aJ and LSV04bl . 

2. Known Results 

2.L The Lattice Packing Problem. The lattice packing problem, arising from the study 
of positive definite quadratic forms, is the oldest and most popular of the three problems 
and has been considered by many authors in the past. As shown in Table [T] the solution to 
the problem was known for dimension d < 8 since 1934. For a description of the extremal 
root lattices Ad, and and the history of the problem we refer the interested reader 
to the book MCS88bl . Recently, Cohn and Kumar IICK04II showed that the Leech lattice A 
gives the unique densest lattice packing in M^'*. Furthermore they showed: The density of 
any sphere packing (without restriction to lattices) in R^^ cannot exceed the one given by 
the Leech lattice by a factor of more than 1 + L65 • 10"'^'^. 



d 


lattice 


density Sd 


author(s) 


2 


A2 


0.9069... 


Lagrange |Lag73|, 1773 


3 


A3 = D3 


0.7404 . . . 


GauB l:Gau40|, 1840 


4 


D4 


0.6168... 


Korkine.Zolotareff I.KZ73I. 1873 


5 


D5 


0.4652 . . . 


Korkine, Zolotareff ||KZ771,1877 


6 


Ee 


0.3729... 


Blichfeldt |Bli34l, 1934 


7 


E7 


0.2953... 


Blichfeldt [Bli34|, 1934 


8 


Eg 


0.2536... 


Blichfeldt IIBli34ll. 1934 


24 


A 


0.0019... 


Cohn, Kumar ||CK04|. 2004 



Table 1 . Optimal lattice packings. 



2.2. The Lattice Covering Problem. The lattice covering problem has only been solved 
up to dimension 5. Recently, we were able to verify the list of known results computation- 
ally. Even more, we found the complete list of 222 local covering optima in dimension 5. 
Table |2] invites to a question formulated by Ryshkov fRys6Tl, who asked for the lowest 
dimension in which A"^ gives not the thinnest lattice covering. In Section|5]we describe the 
method used to find a lattice Lg C M.^ with e(Lg) = 2.4648 . . . < 2.5511 . . . = e(A^). 
Thus the answer to Ryshkov's question is 6. It remains an open problem to prove that the 
lattice L§ gives the best lattice covering in dimension 6. We do not even know the exact 
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coordinates of the lattice yet. Currently, a complete solution in dimension d > 6 seems out 
of reach without completely new methods. 



d 


lattice 


density 0^ 


author(s) 


2 


a; 


1.2091... 


Kershner|Ker39|, 1939 


3 


A3 


1.4635... 


Bambah |Bam54|, 1954 


4 


Aa 


1.7655... 


Delone, Ryshkov |DR63J, 1963 


5 


A^ 


2.1242... 


Rvshkov, Baranovskii iRBTSl, 1975 



Table 2. Optimal lattice coverings. 



2.3. The Lattice Packing-Covering Problem. As in the case of the lattice covering prob- 
lem the lattice packing-covering problem has been solved only for dimensions d < 5. 
And, as in the covering case, we recently were able to verify these results computation- 
ally. Moreover, we found a new best known lattice Lg^ in dimension 6, having a slightly 
lower packing-covering constant 7(_Lg'^) = 1.4110 . . . than the previously best known one 
7(Ea) = V2. 

One reason for studying the lattice packing-covering problem is the open question whether 
there exists a dimension d with 7^ > 2. If so, then any d-dimensional lattice packing with 
spheres would leave space large enough for spheres of the same radius. This would in 
particular prove that densest sphere packings in dimension d are non-lattice packings. This 
phenomenon is likely to be true for large dimensions, but has not been verified for any d so 
far 



d 


lattice 


density 7^ 


author(s) 


2 


A* 


1.1547... 


Ryshkov |Rys74|, 1974 


3 


A^ 


1.2909... 


Ryshkov |Rys74|, 1974 


4 


H4 


1.3625... 


Horvath [Hor82|, 1980 


5 


H5 


1.4494... 


Horvath 0Hor86l. 1986 



Table 3. Optimal lattice packing-coverings. 



3. Lattices and Positive Quadratic Forms 

It is sometimes convenient to switch from the language of lattices to the language of 
positive definite quadratic forms (PQFs from now on). In this section we give a dictionary. 
For further reading we refer to 0CS88bl and IISV04aL 

Given a d-dimensional lattice L = AX^ with basis A we associate a d-dimensional PQF 
Q\x\ — x^A^Ax = x^Gx, where the Gram matrix G = A* A is symmetric and positive 
definite. We will carelessly identify quadratic forms with symmetric matrices by saying 
Q = G and Q[x] = x*Qx. The set of quadratic forms is a (''^^) -dimensional real vector 
space S'^, in which the set of PQFs forms an open, convex cone S^q. The PQF Q depends 
on the chosen basis Aof L. For two arbitrary bases A and B of L there exists aU G G (Z) 
with A = BU. Thus, GLd(Z) acts on S^^ hy Q ^ U^QU. A PQF Q can be associated 
to different lattices L — AX^ and L' = AU^ . In this case there exists an orthogonal 
transformation O with A = OA. Note that the packing and covering density, as well as 
the packing-covering constant, are invariant with respect to orthogonal transformations. 

The determinant (or discriminant) of a PQF Q is defined by dei{Q). The homogeneous 
minimum X{Q) and the inhomogeneous minimum fi{Q) are given by 

X{Q) — min Q[v], = max min Qfa; — v]. 
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If Q is associated to L, then det(L) = ^dct(Q), fi{L) = ^/n{Q), A(L) = y^\{Q)/2. 
Using this dictionary we define 



'^('^^ ^ V ■ 2^' ^ V ■ ^^^^ ^ 1 

We say that a lattice L with associated PQF Q gives a locally optimal lattice packing, 
locally optimal lattice covering or locally optimal lattice packing-covering, if there is a 
neighborhood of Q in 5J^o. so that we have 6{Q) > 5{Q'), <d{Q) < e(Q') or 7(Q) < 
"fiQ') respectively, for all Q' in this neighborhood. 

4. On Packing Lattices 

A PQF Q attaining a local maximum of S{Q) is called extreme. A PQF attaining (5^ = 
maxQg^d^ <^(Q) is called absolutely extreme or critical. One can characterize an extreme 
PQF using the geometry of its minimal vectors 

Min(g) ^{veZ'': Q[v] = A(Q)}. 



Before we state the characterization in Theorem 14.11 we give some more definitions. A 
PQF Q' is called perfect if it is uniquely determined by its minimal vectors, i.e. Q' is the 
unique solution of the linear equations Q[v] = X{Q'), v e Mm{Q'). A PQF Q is called 
eutactic if 

Q^^ e relint cone{ : v e Min(Q)} = < ^ A„t;t>* : > 0, v G Min(Q) 

[ueMm(Q) 

The eutaxy and the polyhedral cone cane{vv*^ : v e Min((3)} also plays an important 
role in the lattice packing-covering problem. As a general reference on basic facts about 
polyhedral cones, which are used throughout this article, we refer to the book of Ziegler 
l|Zie97|. 

Theorem 4.1 (Voronoi IIVor07l ). A PQF is extreme if and only if it is perfect and eutactic. 

This provides an easy way for proving that a given PQF is extreme: after finding the 
minimal vectors, one has to solve a system of linear equations to show its perfectness. 
Then, one has to solve a linear programming problem to verify its eutaxy. By scaling we 
can normalize an extreme PQF Q so that A(Q) is rational. Then, since Q is perfect, the 
matrix entries of Q are rational as well. 

It turns out that there exist only finitely many pairwise non-equivalent perfect PQFs 
in 5>Q. We want to describe two classical algorithms to attain all perfect forms of a given 
dimension d. The first one goes back to Minkowski, the second one is due to Voronoi. Here 
we only state definitions and main results. Additionally, we briefly sketch the computations 
which we were able to perform in low dimensions. We compare them with corresponding 
results in the literature. 

For history and further remarks we refer to IIGL87I . rRB791, fvdWSei, fMarOSI and to 
references therein. Chapter §v of Gruber and Lekkerkerker's book |GL87| gives a compre- 
hensive survey about history, results and literature of the reduction theory of PQFs. The 
article ['RB791 introduces to methods for studying the geometry of PQFs and contains many 
proofs. Van der Waerden's paper [vdW 56 1 is a classic resource for Minkowski's approach. 
The recent book IIMar03l of Martinet gives a contemporary view on Voronoi's approach 
and on possible generalizations. 

4. 1 . Minkowski's Approach. 

Definition 4.2. A PQF Q = {qij ) e 5>o is called Minkowski reduced if 

(i) Q[x] > qu whenever gcd(a;i, . . . , Xd) = 1, i = 1, . . . , d, 

(ii) qi,i+i > 0, i = 1, . . . - 1. 
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Every PQF is equivalent to a Minkowski reduced PQF. The following procedure, which 
is nothing but an algorithmic interpretation of the definition, finds a Minkowski reduced 
PQF equivalent to a given PQF Q. Choose a minimal vector Vi of Q. Then, choose among 
all vectors in H^, which can complement V\ to a lattice basis of Z"^, a vector V2 for which 
the value Q[v2\ is minimal. Using this greedy strategy we get a basis V = (vi, . . . , Vd) G 
GLrf(Z). Now we choose signs for Vi so that v\_iQvi > 0. Hence, V^QV is Minkowski 
reduced. 

The set of Minkowski reduced PQFs forms an unbounded cone in 5>q which is defined 
by the linear inequalities (i) and (ii). By Ai we denote the cone defined by the linear 
inequalities (i) and by A4~^ the one which is defined by the linear inequalities (i) and (ii). 
Minkowski [MinOSl showed that A4, and hence A^+, is a polyhedral cone, i.e. that finitely 
many inequalities (i) imply all others. He showed that every extreme PQF is equivalent 
to a PQF lying on a ray (a one-dimensional face) of Ryshkov | Rys70| proved that 

every perfect PQF is equivalent to a PQF lying on a ray of A^+. On the other hand Cohn, 
Lomakina and Ryshkov [iCLR82J found a ray of C iS>q which contains non-perfect 
PQFs. 

Minkowski [Min87 1 gave a list of conditions implying all others in (i) up to dimension 6. 
Tammela I Tam8 1 j enlarged this list to dimension 7. Besides the d — 1 inequalities qn < 
• ■ • < qdd, the linear conditions for Al, rf = 2, . . . , 7, are attained by plugging the values 
X = {xi, . . . , Xd) E Z"* from Table|4]into (i), where the indices ii, . . . ,id run through all 
permutations of {1, . . . , d}. If d < 7 one has to omit the columns d+1, . . . ,7 and the rows 
with more than d non-zero entries. 





dzX22 














1 


1 



















1 


1 


















1 


1 

















1 


1 
















2 


1 
















1 


1 















2 


1 















2 


2 















3 


2 













1 


1 


1 














2 


1 














3 


1 














2 


2 














3 


2 














2 


2 


2 












3 


2 


2 












4 


2 


2 












3 


3 


2 










2 


4 


3 


2 










3 


3 


2 


2 


2 










4 


3 


2 


2 











Table 4. Tammela's list of linear inequalities defining A4 



We checked these conditions for redundancy using the software 1 r s of Avis fAvi04J . 
Our computations show that in row ^, for d = 6, the entries with xi ^ 3 are redundant 
(as already mentioned by Tammela in IITam73ll ). For d = 7 redundant entries are those 
with xi 7^ 0, 3 and xi — Q,X2 3. In row ^ the entries with xi 3 or X2 7^ 3, and 
in row the entries with a; 1 7^ 4 or 2:2 ^ 3 are redundant. The remaining conditions 
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are all non-redundant and define a facet of A^. In Table |5] we list the number of facets 
and rays as far as we were able to compute them with cdd [Fuk03 1 . We hereby confirm 
earlier results by ° Minkowski liMin87j and Barnes and Cohn [BC76 1. We made the data 
available from the arXiv.org e-print archive. To access it, download the source files for 
the paper arXiv:math.MG/0412320, The files mink2 . ine, . . . , mink6 . ine (due to its 
size the file mink7 . ine is only available from the authors) contain the facets of Ai, the 
files mink2 . ext, . . . , mink 6 . ext contain the rays of A4, the files minkp2 . ine, . . . , 
minkp6 . ine contain the facets of A^+, and the files minkp2 . ext, . . . , minkpS . ext 
contain the rays of A^+. For the data format we chose the common convention (Polyhedra 
_ff -format for the * . ine-files and Polyhedra ^-format for the * . ext-files) of the software 
packages cdd and Irs. 

The computational bottlenecks of Minkowski's approach are apparent. It is not easy 
to find a sufficiently small system of linear inequalities defining A4 (or of A^+). Even if 
one has a minimal system of linear inequalities, then computing its rays is a very difficult 
computational problem in higher dimensions. 



d 




# Facets M 


# Rays M 


# Facets M + 


#Rays M+ 


2 


3 


3" 


3" 


3" 


3" 


3 


6 


12" 


19'' 


9" 


11" 


4 


10 


39" 


323'' 


26'' 


109'' 


5 


15 


200 


15971 


117 


4105 


6 


21 


1675 


7 


1086 


7 


7 


28 


65684 


7 


7 


7 



Table 5. Known number of facets and rays of M and M^. 



4.2. Voronoi's Approach. Now we describe Voronoi's algorithm IIVor07l for finding all 
perfect forms of a given dimension. Let m be a positive number. In the remaining of this 
section we assume that every perfect form Q is scaled so that \{Q) ~ m. The set 

V,n = {Qe S^a : Q[v] > m for all veZ'^\{0}} 

is a convex, locally finite polyhedral cone. Its boundary consists of the PQFs with homo- 
geneous minimum m. A PQF Q is perfect if and only if it is a vertex of Vm- The set of 
perfect PQFs of a given dimension d naturally carries a graph structure which we denote 
as the Voronoi graph in dimension d: Two perfect PQFs Q, Q' are connected by an edge if 
the line segment conv{Q, Q'} (convex hull of Q and Q') is an edge of Vm- In this case we 
say that Q and Q' are Voronoi neighbors. The group GLd{1) acts on Vm, on its vertices 
and on its edges hy Q ^ U^QU . Therefore, one can enumerate perfect PQFs by a graph 
traversal algorithm which we shall now describe. 

Voronoi's first perfect form Q[x\ — X]i=i + '^i<j which is associated to 

the root lattice A^, can serve as a starting point in any dimension. For implementing a 
graph traversal algorithm one has to find the Voronoi neighbors of a given perfect form Q. 
Consider the unbounded polyhedral cone 

V{Q) - {Q' e S'^ : Q'[v] >m,ve Min(Q)}. 

We compute the rays Q + M>oi?i, i = 1, . . . , n of V{Q). The Ri turn out to be indefinite 
quadratic forms. So there are v € Z"^ with Ri\v\ < 0. Then, the Voronoi neighbors of 
Q are Q + piRi where pi is the smallest positive number so that \{Q + piRi) = m and 
Min((5 + PiRi) % Min((5). It is possible to determine pi, for example with the following 
procedure: 
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-pRi 

a+0 
2 



m do 



(a,(3)^ (0,1) 

while Q + PR,^ or X{Q - 

if Q + PRi ^ S^o then /3 ^ 

else (a,/3) ^ (/3,2/3) 

end if 
end while 

while Min(Q + aRt) C Min Q do 

I 2 

if X{Q + -yRi) < m then /3 ^ 7 
else a ^ 7 
end if 
end while 

Pi <^ a 

Perfect forms were classified up to dimension 7. A list of all these forms is given in the 
paper rCSSSa] of Conway and Sloane. One can find an electronic version in the Catalogue 
of Lattice^ by Nebe and Sloane. On his homepag^ , Martinet reports that up to now, 
10916 pairwise inequivalent perfect forms are known in dimension 8 and lists them. 

We verified the results for dimensions < 6 using the programs Irs by Avis |Avi041, 
isomby Plesken and Souvignier I.PS97 1 and shvec by Vallentin [Val99j . In Table|6]we 
give the known classifications of perfect forms, extreme forms and absolute extreme forms 
together with the references where the classifications were established. 



d 




# perfect 


# extreme 


# critical 


2 


3 


1° 


1" 


1" 


3 


6 


1^ 


1" 


1" 


4 


10 


2" 


2= 




5 


15 


3"* 


3-^ 


1^ 


6 


21 


If 


6"= 


1/ 


7 


28 


33'' 


30'' 




8 


36 


> 10916^ 


? 




9 


45 


> 500,000 


? 


? 


24 


300 


7 


7 


V 



Lagrange |Lag73 1 

GauB IIGau40l 

Korkine, Zolotareff fKZ73l 
Korkine, Zolotare ff LKZTTl 
Hofreit er |lHo0 3| 
Barnes rBar57| 
Vetchinkin I Vet82l 



Jaquet-Chiffelle llJC93l 

Laihem, Baril, Napias, Batut, Martinet (see http : // w w w.math .u-bordeaux. fr/~martinel) l 
Cohn, Kumar|CK04J 



Table 6. Classifications of perfect, extreme and absolute extreme forms 



The computational bottleneck of Voronoi's approach is mainly the enumeration of all 
rays of the polyhedral cone V{Q) in case of a large set Min(Q) of minimal vectors. Mar- 
tinet writes in ||Mar031 . Ch. 7.1 1 : "The existence of Eg [...] makes hopeless any attempt to 

'|http://www.research.att.comr njasAattices/pe rfect.html | 
^ [http://www.math. u-bordeaux. t'r/~martinet^ 
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construct the Voronoi graph in dimension 8". Another problem is the combinatorial explo- 
sion, when d > 9. We found more than 500000 inequivalent perfect forms in dimension 9 
and we strongly believe there exist millions of them. 

Finally, we want to remind of Coxeter's A^-hypothesis: Although finding perfect forms 
with maximal packing density is a very difficult problem, finding perfect forms with min- 
imal packing density might be very easy. In IICox51l Coxeter formulates the following 
conjecture: 

Conjecture 4.3. (Coxeter's A^-hypothesis) Voronoi's first perfect form gives the minimal 
packing density among all perfect forms of a given dimension. 

Our computations support Coxeter's conjecture. But on the contrary, Conway and 
Sloane |.CS88aJ conjecture that it is false for sufficiently large d. 

5. Covering and Packing-Covering Lattices 

The lattice covering problem and the lattice packing-covering problem can be treated in 
parallel. We describe below that both problems have only finitely many local optima which 
can be found by solving finitely many convex optimization problems. This is mainly due 
to Voronoi's theory of Delone subdivisons, which we briefly review. For a detailed account 
we refer to [SV04ail . 

With an implementation of the proposed algorithms we found all local optima in dimen- 
sion d < 5 and some new best known lattices in dimension d > 6. 

For both problems, recognition of local optima is not as easy as for the lattice packing 
problem. Due to the involved convexity we can give sufficient conditions for local optima, 
allowing to compute a certificate for the local optimality of a lattice. This is in particular 
applicable, if the Delone subdivison is a triangulation which is the generic case (for defini- 
tions see below). Exemplarily we give a proof of the local packing-covering optimality of 
the lattices AJ. 

In some cases it is possible to attain good or even tight "local lower bounds" for the 
lattice covering density and the packing-covering constant. This is demonstrated for the 
local packing-covering optimality of the Leech lattice. A similar proof of the local covering 
optimality of the Leech lattice is given in [SV04b|. 

5.L Voronoi's Theory of Delone Subdivisions. Let Q be a positive semidefinite qua- 
dratic form. A polyhedron P = conv{t)i, t)2, . . .} with Vi,V2, ■ ■ ■ G TL'^, is called a DeZone 
polyhedron of Q if there exists a c e M'' and a real number r e M with Q\vi — c\— for 
all z = 1, 2, . . ., and Q\v - c] > for afl other v e Z"* \ {vi, V2, • ■ •}■ The set Del(Q) 
of all Delone polyhedra is called the Delone subdivision of Q. It is a periodic face-to-face 
tiling of M"^. Therefore Dc\{Q) is completely determined by all Delone polytopes having 
a vertex at the origin 0. We call two Delone polyhedra L, L' equivalent if there exists a 
V G Z'^ so that L = V ± L' . Note moreover that the inhomogeneous minimum is 
at the same time the maximum squared circumradius of its Delone polyhedra. We say that 
the Delone subdivision of a positive semidefinite quadratic form Q' is a refinement of the 
Delone subdivision of Q, if every Delone polytope of Q' is contained in a Delone polytope 
ofQ. 

By a theory of Voronoi llVorOSII . the set of positive semidefinite quadratic forms with a 
fixed Delone subdivision D is an open (with respect to its affine hull) polyhedral cone in 
5>Q. We refer to this set as the secondary cone A(I?) of the subdivision. In the literature 
the secondary cone is sometimes called L-type domain of the subdivision. The topological 
closure A (I?) of a secondary cone is a closed polyhedral cone. The relative interior of each 
face in S^q is the secondary cone of another Delone subdivision. If a face is contained in 
the boundary of a second face, then the corresponding Delone subdivision of the first is a 
true refinement of the second one. 

The interior of faces of maximal dimension C^^^) contain PQFs whose Delone subdivi- 
sion is a triangulation, that is, it consists of simplices only. We refer to such a subdivision 
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as a simplicial Delone subdivision or Delone triangulation. As mentioned in Section [3] 
the group GLrf(Z) acts on 5>q. One of the key observations of Voronoi is that under this 
group action there exist only finitely many inequivalent Delone subdivisions, respectively 
secondary cones. 

Theorem 5.1 (Voronoi BVorOSII '). The topological closures of secondary cones of Delone 
triangulations give a face-to-face tiling of 5>q. The group GLrf(Z) acts on the tiling, and 
under this group action there are only finitely many non-equivalent secondary cones. 



(1) 



minimize 
subject to 



Given a Delone triangulation V, the Delone triangulations V with A (I?') sharing a 
facet with A(I?) are attained by bistellar operations (flips). These change a triangulation 
only in certain repartitioning polytopes associated to the facet. By this operation it becomes 
possible to enumerate all Delone triangulations, and hence all Delone subdivisions in a 
given dimension. For details we refer to IISV04al . 

5.2. Obtaining Local Optima via Convex Optimization. For a fixed triangulation P, we 
can formulate the lattice covering, as well as the lattice packing-covering problem in the 
framework of Determinant Maximization Problems. Following Vandenberghe, Boyd, and 
Wu IIVBW98I their general form is 

c*a; — logdet G{x) 
G{x) ^ 0,F{x) h 0. 

Here, the optimization vector is a; e M^. The objective function contains a linear part 
given by c e and G : ^ M™^", F : ^ R"^" are both affine maps 

G{x) = Go+xiGiH \-xdGd, 

F{x) = Fo+xiFi + --- + xdFd, 

where Gi G K^x"*^ i?. ^ R"X", i = 0, . . . , Z?, are symmetric matrices. The notation 
G{x) >- and F{x) y gives the constraints "G(a;) is positive definite" and "F{x) is 
positive semidefinite". Note that we are dealing with a so-called semidefinite programming 
problem, if G{x) is the identity matrix for all x G MP . 

For the lattice covering, as well as the lattice packing-covering problem, we can express 
A'(Q) < 1 as a linear matrix inequality F{qij) y with optimization vector Q — {qij). To 
see this, it is crucial to observe that an inner product (•, •) defined by {x, y) — x*Qy gives 
a linear expression in the parameters (q^ ) for any fixed choice a;, y e Z'' (or R''). Delone, 
Dolbilin, Ryshkov and Stogrin ODDRS70I showed 

Proposition 5.2. Let L = conv{0, Vi, . . . , Vd} C R'' be a cZ-dimensional simplex. Then 
L's circumradius is at most 1 with respect to (•, •) if and only if 

/ 4 {vi,Vi) {V2,V2) ... {Vd,Vd)\ 

{vi,Vi) {vi,Vi) {vi,V2) ... {vi,Vd) 
{V2,V2) {V2,vi) {V2,V2) ... {v2,Vd) 



BRl(Q) 



\iVd,Vd) {Vd,Vi) {Vd,V2) 



>- 0. 



iVd,Vd)J 



Since a block matrix is semidefinite if and only if the blocks are semidefinite, we con- 
clude 

Proposition 5.3. Let Q — (qij) e S'^q be a PQF Let I? be a Delone triangulation refining 
Dc\{Q), and let Li, . . . , L„ be a representative system of all non-equivalent d-dimensional 
Delone polytopes in T). Then 

/|BRr„(Q) I ... \ 



KQ) < 1 



BRli(Q) 






BRl,(Q) 






BRlJQ) 



0. 
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Thus < 1 can be brought into one linear matrix inequality of type F{qij) >z 0. 

We can moreover add linear constraints on the parameters qij by extending by a 1 x 1 
block matrix for each linear inequality. In this way we can get one linear matrix inequality 
for the two constraints ij,{Q) < 1 and Q £ A(I?). 

For a fixed Delone triangulation V, we can therefore determine the optimal solutions of 
the lattice covering problem of all PQFs Q for which 2? is a refinement of Del(Q). Recall 

that the covering density of a PQF Q in d variables is G(Q) = \J jc'^q ' '^d- Scaling of 
Q by a positive real number a leaves 8 invariant. Thus we may maximize det{Q) while 
f^iQ) < 1- For all Q £ A(I?) this can be achieved by solving 

minimize — log det (Q) 
subject to Q :^ 0, 

Q £ AiV), /i(Q) < 1. 

With an analogues specialization of problem ([T]), we are able to attain optimal solutions 
of the lattice packing-covering problem among all PQFs Q for which 2? is a refinement 
of Del((5). Because 7(Q) = 2 • VM^QV^MQ)' '^^ have to maximize X{Q) while < 
1. Maximizing X{Q) is not as straightforward as maximizing the determinant, since we do 
not know which vector is the shortest. By a theorem of Voronoi llVorOSI we know though 
that among the (at most 2(2"^ — 1)) edges [0, v] = conv{0, v} £ D, there exist some with 
V £ Mm{Q). Consequently, if we require Q[v] < Q[w] for all w with [0,w] £ T), we 
know V £ Mm{Q), respectively X{Q) = Q[v]. So we can solve for each [0, v] £ T) the 
semidefinite programming problem 

minimize —Q[v] 

subject to Q[v] < Q[w], for all w with [0, w] £ V, 

g e a(p),m(Q) < 1- 

Note that in many cases the constraints have no feasible solutions, since in general not 
all of the V with [0,v] G I? are elements of Min(Q). 

So both, the lattice covering as well as the lattice packing-covering problem, are reduced 
to convex programming problems if restricted to the closure of a secondary cone. Conse- 
quently, there is at most one local minimum of O, respectively 7, for each of the cones 
A(I?). This was first observed by Barnes and Dickson IIBD67II in the covering case and by 
Ryshkov ]Rys74[ for the packing-covering problem. 

Proposition 5.4. Let P be a Delone triangulation. Then there exists a unique minimum of 



(1) {0(Q) : Q £ A(I?)} and the set of PQFs attaining it is equal to all positive 
multiples of a single PQF. 

(2) {^{Q) : Q £ A(P)} and the set of PQFs attaining it is convex. 

In case of a triangulation T), it follows that Q £ A (I?) is a locally optimal solution with 
respect to 8 or 7 if and only if it is an optimal solution within A(I?). In general we have 
the trivial 

Proposition 5.5. A PQF Q is a locally optimal solution with respect to 9 or 7, if and only 
if it is an optimal solution for all Delone triangulations T> with Q £ A(I?). 

5.3. Computational Results. By the foregoing propositions we know that the number of 
local optima is bounded from above by the number of pairwise inequivalent Delone tri- 
angulations in M.'^. Voronoi MVorOSII classified these triangulation in dimension 2, 3 (only 
one each) and 4 (three). By the work of Baranovskii and Ryshkov IIBR73L MRB78I En- 



gel |Eng98|, and Engel and Grishukhin |EG02j we know of exactly 222 Delone trian- 
gulations in dimension d = 5. Using Irs BAvi04l and an implementation (in C++) of 
Voronoi's algorithm for enumerating Delone triangulations, we were able to confirm these 
results tSV04aJ . For dimension d > 6 we experience a combinatorial explosion, e.g. Engel 
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I Eng04| reports on more than 2, 129, 120 pairwise inequivalent Delone triangulations for 

d = 6. 



d 




# covering optima 


# packing-covering optima 


2 


3 


1" 


1'' 


3 


6 




1" 


4 


10 


3^^ 


3'= 


5 


15 


222^ 


47/ 



° Kershner IKer39ll 
^ Ryshkov |jRys74| 
" Bambah |Bam54| 



Baranovski [ Bar65 L Dickson IIDic67l 
" Horvath IHor82l 
f Schiirmann, Vallentin OSV04al 

Table 7. Classifications of locally optimal covering and packing- 
covering lattices. 



For each of the triangulations in dimension d < 5 we determined the local optima with 
respect to 8 and 7 using the software package MAXDET^of Wu, Vandenberghe, and Boyd 
as a subroutine. By this we confirmed the known results for dimensions d = 2,3,4 and 
extended them to d = 5 (see Table|7ll. Note that for the lattice covering problem there exists 
a local optimum for each triangulation, while this is not the case in higher dimensions (see 
IISV04bll ) and for the lattice packing-covering problem. 

Using our implementation we also found two lattices which currently give the best 
known covering and packing-covering in dimension 6: 

Theorem 5.6 ( IISV04al ). In dimension 6, there exits a lattice Lg with 6(L§) = 2.4648 . . . 
and a lattice with -/{LP") = 1.4110 . . .. 

In IISV04bll we show that the root lattice Eg does not give a locally optimal lattice cov- 
ering, by constructing a refining triangulation T) of Dc1(Qe8 ) in which 8's local optimum 
is not attained by the PQF Qes- The PQF found in this way even beats the formerly best 
known value 0(Ag) by more than 12%. By looking at a bistellar neighbor of the triangula- 
tion V, we found the currently best known covering lattice in dimension 8. 

Theorem 5.7 ( ||SV04bl ). In dimension 8, there exists a lattice Lg with9(Lg) 3.1423 . . .. 

Looking at the results in dimension d = 6, 8 it is interesting to observe that we found the 
new covering lattices by looking at triangulations refining the Delone subdivisions of the 
lattices E^. By looking at a corresponding refinement of E^, we also found a new covering 
record in dimension 7. It remains to see if these results have a common explanation. . . 

5.4. Sufficient Conditions for Local Optima. A disadvantage of finding local optima via 
convex programming is that solutions can only be approximated. But this is an inherent 
problem: In contrast to the lattice packing problem, local optima to the other two problems 
can in general not be represented by rational numbers. One has to use algebraic num- 
bers instead. In some cases it might be possible, e.g. with additional information on the 
automorphism group, to attain exact coordinates from a first approximation (see IISV04all 
for an example). In other cases we might have a conjectured optimal form Q' and want 
to compute a "certificate" verifying its local optimality. The following two propositions 
give such a criterion in terms of the gradient ghiQ') = grad | BR^ \ {Q') of the regular 
surfaces | BRl((3)| = at Q' . Both are a consequence of the geometric fact that we 

- ^http://www.stanford.edurboyd/MAXDET.fi5iI] 
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have a local optimum at Q' e A(P) with fJ,{Q') 1 if and only if there exists a hyper- 
plane through Q', separating the convex sets {Q G S'^q : L G | BRl((5)| > 0} and 
{Q £ : det(Q) > det{Q')}, respectively {Q G : A(Q) > X{Q')}. 

Proposition 5.8 (Barnes and Dickson IIBD67I ). Let V he a Delone triangulation. Then 
Q G A (I?) with fi{Q) = 1 is a unique locally optimal solution to the lattice covering 
problem if and only if 



Q- 



'Cone{5L(Q) : L G P with | BRl(Q)| = 0}. 



The corresponding result for the lattice packing-covering problem is 

Proposition 5.9 ( IISV04al ). Let I? be a Delone triangulation. Then Q G A(I?) with 
= 1 is a locally optimal solution to the lattice packing-covering problem if and only 

if 

cone{vv* : v € Mm{Q)} n - cone{gL{Q) : L G P with | BRl(Q)| = 0} 7^ 0. 

Combining these two propositions we get 

Corollary 5.10. Let I? be a Delone triangulation. Then Q G A (I?) is a unique locally 
optimal solution to the lattice packing-covering problem, if Q is eutactic and a locally 
optimal solution to the lattice covering problem. 

Example 5.11. We can use Corollary EJO] to show that A^, d > 2, with j{A*j) = 
is locally optimal for the lattice packing-covering problem. Ryshkov | |Rys74[ gave another 
proof of this fact. The lattice is known to give a locally optimal lattice covering (see 
IIGam62ll . IIGam63l . IIBle62ll ). A quadratic form Qa* associated with AjJ is 



with Ai(QA*) 



d{d+2) 
12 



and A(g) 



A 
d. 



-1 



-1 



-1\ 



-1 
d/ 



The set Min((5Aj) contains exactly 2{d 



elements, namely the standard basis vectors ei, . . . . e^, their negatives and ± J2i 
(see IICS88bll ). Thus in particular 



1) 

e,: 



{d+l)-Q^} = Q^, = 



/2 
1 

Vi 



1 



1\ 



1 

2/ 



E 

T)T)':i;eMin(QA« ) 



VV 



Hence, Qa* is eutactic because Q/^} G relint cone{t>t>* : v G Min(QA*)} and therefore 
the assertion follows. 



Propositions 15.81 and 15.91 assume that V = Del(Q) is a Delone triangulation. If this is 
not the case, the situation becomes more complicated, in particular for the lattice packing- 
covering problem. 

For the covering problem we only have to add a condition on the set 

IJ {Q G A{V') : |BRl(Q)| > OforallL G V'}. 

'D'<T> 

where V < V denotes that V is a Delone triangulations refining V. This set is a subset 
of {Q G : /i(Q) < 1}. We require that is sepamtable at Q, that is, there exists 
a supporting hyperplane of Vt> through Q. This is in particular the case, if there exists a 

small r > such that {Q + rB^ ^ ') r\ Vv is convex. 

Proposition 5.12. Let 2? be a Delone subdivision and Q G A(I?) with ^(Q) = 1. Then 
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(1) Q is a locally optimal solution to the lattice covering problem, if and only if Vx> is 
separatable at Q and 

Q-^ e - cone{5L(Q) : i G 2?' < 2? with | BRl(Q)| = 0}. 

(2) Q is a locally optimal solution to the lattice packing-covering problem, if Vv is 
separatable at Q and 

cone{vv* : v e Min(Q)} n - cone{gL(g) : LeV < P with |BRl(Q)| = 0} ^ 0. 

In case of the lattice packing-covering problem the "only if" part is missing, because we 
can not exclude the case of a locally optimal solution Q'with not being separatable at 
Q'. This is due to the fact that {Q E S^q : X{Q) > X{Q')} is not smooth in contrast to 

{Q G : det(Q) > det(Q')}- 

This phenomenon seems to happen to PQFs associated to the root lattice Eg. This lattice 
is known to give a globally optimal solution to the lattice packing problem, but not a locally 
optimal solution to the lattice covering problem (see llSV04bl ). Nevertheless computational 
experiments support the 

Conjecture 5.13. The root lattice Eg gives a locally optimal solution for the lattice packing- 
covering problem. 

Zong OZon02l even conjectured that Eg gives the unique globally optimal solution to the 
lattice packing-covering problem in dimension 8. 

5.5. Local Optima via Local Lower Bounds. In OSV04al we describe a way to attain 
local lower bounds for the covering density and the packing-covering constant due to 
Ryshkov and Delone. A variant of this method is successfully used in [SV04bl to prove 
the local covering optimality of the Leech lattice. Here we describe a corresponding local 
lower bound for the lattice packing-covering problem. As an example we use it to prove 
the local packing-covering optimality of the Leech lattice directly. 



Proposition 5.14. Let Li , . . . , L„ be a collection of Delone simpUces of a PQF Q. Then 

7(Q)>2^ 



/ trace(i^ • Qp) 



{d+l)X{QF) 

with the PQF F = J2^ Ek^i ^^.kvll and a PQF Qf with 

F e conc{t;t;* : v e Mm{Qp)}. 

One can prove this Proposition by doing obvious modifications to the proof of Propo- 
sition 10.6 in MSV04al . As in Proposition 15 .9 1 we use the following fact: A linear function 
f{Q) = trace(F • Q), with a PQF F, has a minimum on the homogeneous minimum A 
surface {Q G S^q : X{Q) — A} at Qp if and only if e cone{vv* : v E Mm{Qp)}. In 
particular, if Q is eutactic and F = Q^^, then Proposition 15 . 1 41 is immediately applicable 
with Qp = Q. 

Example 5.15. We use Proposition l5. 14l to show that the Leech lattice is a locally optimal 
packing-covering lattice. 

Let us briefly review some necessary properties of the Leech lattice A. For further 
reading we refer to |CS88F|. An associated PQF Qa has (up to congruences) 23 different 
Delone polytopes attaining the maximum squared circumradius /^(Qa) = 2. One of them 
is the Delone simplex L of type A24. 

Now we apply Proposition 15 . 1 41 to the orbit of L under the automorphism group Coo = 
{T e GL24(Z) : T*QaT = Qa} of Qa. We get F = EtsCoo Ee ee*, where e 

runs through all the edge vectors of TL. In ||SV04bl it was shown that F = Qa^- 
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to the fact that Min Q\ is a spherical 2-design with respect to the inner product given by 
Qa, we know (see BSV04bll for details) 

DeMin(QA) 

Thus, Qa is eutactic and we may use Qf = Qa in Proposition 15. 141 With A(Qa) = 4 we 
derive 

OL . 24 

for all PQFs Q with Delone simplices TL, T e Coq, which proofs the assertion. 
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